Optical Precursors with Self-induced Transparency 
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Optical Sommerfeld-Brillouin precursors significantly ahead of a main field of comparable am- 
plitude have been recently observed in an opaque medium with an electromagnetically induced 
transparency window [Wei et al, Phys. Rev. Lett. 103, 093602 (2009)]. We theoretically ana- 
lyze in this paper the somewhat similar results obtained when the transparency is induced by the 
propagating field itself and we establish an approximate analytic expression of the time-delay of the 
main-field arrival, which fits fairly well the result obtained by numerically solving the Maxwell-Bloch 
equations. 

PACS numbers: 42.25.Bs, 42.50.Md, 42.50.Gy 



More than one century ago, Sommerfeld examined the 
apparent inconsistency between the existence of superlu- 
minal group velocities and the theory of relativity. Con- 
sidering an incident field switched on at time t — (step 
pulse), he showed that, no matter the value of the group 
velocity, no field can be transmitted by a linear disper- 
sive medium before the instant t = L/c where L is the 
medium thickness and c the velocity of light in vacuum 
[l|. Subsequently he and Brillouin studied the fast os- 
cillatory transients appearing at t > L/c in the partic- 
ular case of a single-resonance Lorentz medium 0, 0]. 
They named them forerunners insofar as, in proper con- 
ditions, they can distinctly precede the establishment of 
the steady-state field (the main field). Renamed optical 
precursors, forerunners have entered classical textbooks 
[4|, |5J] and continue to raise a considerable interest. The 
theoretical results of Sommerfeld and Brillouin have been 
improved, even rectified (the amplitude of the precursors 
was in particular strongly underestimated in their work) , 
and different models of linear dispersive media have been 
considered. See for a recent review. 

Despite the abundant literature on precursors, there 
are very few papers reporting direct observation of pre- 
cursors distinguishable from the main field. The diffi- 
culty of such an observation has been soundly discussed 
by Aavikssoo et al. [3] who achieved in 1991 an exper- 
iment involving single-sided exponential pulses (instead 
of step-pulses) and exploiting the dispersion originating 
from a narrow exciton line in AsGa Q . For proper detun- 
ing of the optical carrier frequency uj c from the resonance 
frequency ujo, optical precursors appear as a small spike 
superimposed on the main pulse. See also The ob- 

servation of precursors significantly ahead of a main field 
of comparable amplitude obviously requires the use of 
long enough square pulses and of a medium fairly trans- 
parent at the optical carrier frequency, the corresponding 
group delay being long compared to the duration of the 



precursors. As discussed in |lll.ll2l|. the latter conditions 
are met in an opaque medium with a narrow transparency 
window (slow light medium). Such an experiment has 
been recently achieved by Wei et al. [l3| in an opaque 
cloud of cold atoms with an electromagnetically induced 
transparency (EIT) window. Note that, in this experi- 
ment (as in all the studies of precursors), the propagating 
field linearly interacts with the medium. 

For comparison, we will examine here the nonlinear sit- 
uation where the medium transparency is induced by the 
propagating field itself (13 . [H| • Figure [TJ shows the re- 
sult of an experiment achieved in such conditions [l6j |. 
The medium is a gas of HC 15 N at low pressure con- 
tained in a 182m-long oversized waveguide and the inci- 
dent wave is on resonance with the molecular rotational 
line J = 0, M = -> J = 1, M = (wavelength 
A c « 3.5 mm). The gas behaves as a 2-level medium 
[r7| characterized by T\ (T2 ) the relaxation time for the 
population difference (the polarization), T 2 * the Doppler 
time and a the resonant absorption coefficient at low in- 
tensity (extrapolated from the Lorentzian wings of the 
line). See [18] for details. The incident wave is charac- 
terized by Iq its intensity normalized to the saturation in- 
tensity and r r its rise time. The observed step responses 
clearly have some similarities with those obtained in the 
EIT experiment with a short transient preceding the 
establishment of a steady state regime (main field) . The 
quasi Rabi oscillations [15| accompanying the latter are 
obviously absent in the EIT experiments but oscillations 
having a linear origin (postcursors) can also be observed 
in this case [12]. 

To analyze the previous results, we provisionally ne- 
glect the Doppler broadening and assimilate the guided 
wave to a plane wave propagating in the z direction 
(0 < z < L), with an electric field polarized in the x 
direction. As long Ti, T 2 » 1/oJc and a <C w c /c, 

the sl owly varying envelope approximation (SVEA) [l7j 
holds (l9| and we write the E x component of the electric 
field as : 
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E x (z,t)=Re 
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Figure 1: Observed step response of a resonant absorbing 
medium. Parameters : aL ~ 200 , Ti ~ Ti « 10 /xs {TijaL ps 
50 ns), T 2 * » 1.3 /xs , w c = 5.4 x lO 1 ^ 1 (w" 1 « 1.8 ps ), 
r r = 12ns ; J ~ (a) 2100,(b) 960 (c) 620 (d) 350. In each case, 
the intensity is normalized to that of the step transmitted in 
the absence of gas (dashed line). 



where, as in all the following, t is a local time (real time 
minus z/c), and E(z,t) is the slowly varying field en- 
velope. Denoting u the dipole matrix element for the 
transition (chosen real), R(z,t) — fj,E(z,t)/h the Rabi 
frequency, n(z, t) the population difference per volume 
unit (no its value at equilibrium) and P(z, t) the envelope 
of the electric polarization induced in the medium, it is 
convenient to introduce the dimensionless quantities D = 

n /"o, P = ^7^5 and E = ^E^fTm/h = Ry/T^, 

all real in the resonant case. I = E 2 is the intensity nor- 
malized to the saturation intensity. The Maxwell-Bloch 
(MB) equations governing the system evolution take then 
the simple form 

9E (2) 



dz 



-P 



dP 

T 2 —=DE-P 
at 



(3) 



response E(L,t) of the medium (with a time resolution 
equal to r r ) is then obtained by solving the MB equations 
with P(z,0) = 0, D{z,0) = 1 and E(0,t) = E Q(t) = 
\/7o0(f) where 0(t) is the unit step function. This prob- 
lem has been examined by Crisp [15| when the relaxation 
effects are negligible, a condition obviously not met in the 
experiments. 

The long term behavior of the step response (t ^ 
T\,T 2 ) is obtained by solving the MB equations in steady 
state. Combining EqsHandH, we find P = E/(l + E 2 ) 
and, putting this result in Eq[2l we easily retrieve the 
transmission equation [2p| . l2ll . |22 | 



7(oo) + ln/(oo) = I +1nI - aL 



(5) 



where I(t) is a short hand notation of the transmitted 
intensity I(L,t). The medium being optically thick in 
the linear regime (aL ^> 1), the absorption is fully sat- 
urated (i(oo)//o « 1) only when the incident (normal- 
ized) intensity is extremely large (Iq ^> aL). In fact 
the transmitted field (main field) will be significant (par- 
tial transparency) as soon as lo — aL — O(aL). The 
transmission equation takes then the approximate form 
I(oo)/Iq » 1 — aL/ Iq and a transmission I(oo)/Iq > 1/3 
is obtained for Iq > ZaL/2. 

Consider now the short term behavior of the step re- 
sponse. By combining the integral form of Eqs[3]and[4] 
and taking into account that D(z,t) < 1, one can estab- 
lish the inequality [23] 



1-D(z,t) < 



R(z,t')dt' 



< R 2 t 2 



(6) 



where Rq is the Rabi frequency associated with the in- 
cident step ( Rl = ijhjr). When R^t 2 < 1, D(z,t) ~ 1 
and the MB equations are reduced to the couple of lin- 
ear equations dE/dz = -aP/2 and T 2 dP/dt = E - P. 
So, at least in this time domain and though Iq ^> 1, the 
medium behaves as a linear system (small pulse-area ap- 
proximation [HI])- Its response E(L, t) is easily retrieved 
from the previous couple of equations and can be written 
asfHSl 



E{L,t) = E O{t) 



t/T 2 - 



1-aL 



(V2aLuj 
y2aLu 



e~ u du 



(7) 

When aL 1 , the integral can be transformed to obtain 
E(L,t) » E Q(t)e- t/T2 J (y/2aLt/T 2 ) (8) 



dD 

T 1 — = -PE+(l-D) 
at 



(4) 



We assume that the rise time r r of the incident in- 
tensity, while long compared to l/co c (as above men- 
tioned), is short with respect to all the other character- 
istic times of the system (1/R, Ti, T 2 and T 2 /aL). The 



For x > 1, J (at) ~ y ™ cos { x ~ f) and E(L,t > 
g£) « E+(t) + E-(t) where 



E± (t) = -5= e -*/ T " 



exp 



±i {^/2aLtjT 2 - tt/4 



2tt 



(2aLt/T 2 ) 



i/4 



(9) 
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So the optical field is made of two components of equal 

amplitude and instantaneous frequency uj c ± yj $jr, 

which are nothing but that the Sommerfeld (-E+) and 
Brillouin (E-) precursors as determined by the saddle 
point method of integration [H, [H, ■ The linear char- 
acter of the short-term response (and thus its analysis 
in terms of precursors) is well supported by the exper- 
iments. As shown Fig[H the shape of the correspond- 
ing transient is roughly independent of the incident in- 
tensity. By numerically solving the MB equations, we 
find that the condition R\t 2 <§; 1 is much too severe and 
that the linear approximation satisfactorily holds up to 
t = 2tt/Rq, the Rabi period of the incident field. It even 
holds later in the experiments because the transversal 
inhomogeneity of the field partially washes out the (non- 
linear) quasi Rabi oscillations while it does not affect the 
linear response (the precursors). 

In the EIT experiments, the probe field linearly inter- 
acts with the medium at every time and the arrival of the 
main field is determined by the (slow) group velocity [12]. 
In the present case, this arrival is fixed by fully nonlinear 
phenomena, the study of which requires the resolution 
of the complete MB equations. We first examine the 
solution obtained in the rate equations approximation 
(RE A) The equations to solve are then reduced to 

dl/dz = -alD and T x dD/dt = -U((l + 1) + 1 (23,11 
with D(z,0) = 1 and 1(0, t) = I @(t). Eliminating D 
and integrating in z, we get [22|] 



Tid(lnJ) fdt = lnJ + I - aL - In I - I 



(10) 



with 1(0) = Iq exp (—aL). The transmitted intensity I(t) 
is finally given by the implicit equation 



t 



lit) 



dV 



Iq cxp(— aL) 



F (In 7 + J - aL - In I' - 1') 



(11) 



The transmission T(t) = I(€)/Iq monotonously increases 
from exp (— aL) to I(oo)/Io, where I(oo) is given by 
EqUJ In the conditions considered here [aL ^> 1, 
Io-aL = O(aL)], T(0) « 0, T(oo) w 1 - aL/I Q and 
the transition between these two values is very steep 
(Figl2]). The time-delay of the arrival of the main in- 
tensity is conveniently defined as the time Td such that 
7(r d ) = 7(oo)/2. It is given by EqUUby taking J(oo)/2 
as upper limit of integration. When Iq ^> aL (full sat- 
uration limit), EqQT]can be explicitly integrated to give 
T(t) f» [1 + ex p (a L — I t/Ti)]~ l in agreement with the 
result given in [2a|. The 10-90% rise time At of the inten- 
sity and the time-delay Td then read as At « 4 In 3 (Ti/Iq) 
and Td ~ aL (T\/I$) <C T\. When the saturation is only 
partial, the time-delay Td as a function of increases 
much faster than aLT\ / Iq and values of the order of T\ 
can be attained while keeping a significant transmission 
(Fig©. 

The REA does not take into account the coherent ef- 
fects. It eliminates in particular the quasi Rabi oscilla- 
tions accompanying the main field. One may however 
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Figure 2: Step response obtained in the rate equations ap- 
proximation (REA) as a function of the normalized time t/T\ , 
Optical thickness aL = 200. Each step response is labeled by 
the corresponding incident intensity Jo- The step response 
obtained for Jo —* oo is given for reference (dashed line). 



expect that the signals obtained by this way are a satis- 
factory approximation of the exact signals, the oscillatory 
parts of which would have been filtered out. To check this 
idea, we have compared, aL and Iq being fixed, the step 
response obtained by using the REA (independent of T2) 
to those obtained by numerically solving the MB equa- 
tions for two different values of T%jT\ (Fig|3|). The three 
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Figure 3: Comparison of the step response obtained by the 
REA (dashed line) to those obtained by numerically solv- 
ing the Maxwell-Bloch (MB) equations with T 2 /T l = 1 (full 
line) and with T2/T1 = 1/2 (dotted line). Other parame- 
ters: aL = 200 and I = 400, leading to T(oo) w 1/2. Note 
that the pseudo-period of the oscillations superimposed to the 
steady state in the MB solutions are nea rly equal to the cor- 
responding Rabi period 2tt/Rcc = 2ty^/TiT2 / 'I (00), namely 
0.44Ti (0.31Ti) for T2/T1 = 1 (T2/T1 = 1/2). 

step responses are obviously different but the time-delays 
Td (as defined before) are very close. Similar simulations 
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made for different values of the parameters show that this 
result is not accidental. It appears that EqfTTI provides 
the exact time-delay with a precision better than 10% in 
all the cases of physical interest, that is when the pre- 
cursors are well developed before the arrival of the main 
field and the latter has a significant amplitude. 



1.0- 
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Figure 4: Numerical solution of the MB equations taking into 
account the transverse distribution of the field (full line) for 
aL = 200, Io = 400, and T 2 = Ti. The REA (dashed line) 
and MB (dotted line) solutions obtained with the plane wave 
model are given for reference. 

We will now examine the modifications brought to the 
step response by some effects neglected in the previous 
theoretical analysis. The most important one results 
from the transverse inhomogeneity of the guided wave. 
Figure [4] shows a typical step-response obtained by us- 
ing a MB numerical code extended to include a trans- 
verse variation of the field As expected, the linear 
part of the response (precursors) is not changed (it is 
even slightly prolonged) but the quasi Rabi oscillations 
(strongly depending on the field amplitude) are dramat- 
ically affected. Their amplitude is considerably reduced 
and their damping is accelerated, in agreement with the 
experimental result (FigJT]). However we remark that 



the time-delay tj, is not significantly larger than that ob- 
tained in the plane-wave and rate-equations approxima- 
tions. Similar calculations including the Doppler broad- 
ening instead of the field inhomogeneity in the plane- wave 
MB numerical code show that, even when T 2 * = 0.13T2 
(parameters of FigdJ, the Doppler effect negligibly af- 
fects the precursors and slightly reduces the time-delay 
Trf. This can be explained by observing that the right 
time scale for the precursors and the nonlinear response 
is not T 2 but, respectively, T 2 /aL <C T 2 * and l/i?o < T 2 . 
Finally, the finite rise time of the incident step essen- 
tially affects the most rapidly varying part of the step 
response, namely the transient associated with the pre- 
cursors and first the intensity h of its first peak. When 
aL ^> 1, Ji only depends onr — aLT r /T 2 and attains the 
intensity Iq of the incident wave when r <C 1. This con- 
dition is approximately met in the experiment reported 
in [l6|] where I\ ~ 0.9/o (Fig[T|). Similar results could be 
obtained at optical wavelength by propagating a Gaus- 
sian beam in an ensemble of laser cooled 2-level atoms. 
We have then T 2 ^> T 2 and the Doppler effect negligi- 
bly affects the precursors and the quasi Rabi oscillations. 
In other respects, T 2 (typically 50ns) is about 200 times 
shorter than in the microwave experiment. For a good 
observation of the precursors, the rise time of the inci- 
dent step should also be 200 times shorter, namely in the 
50ps range (attained with electro-optic modulators). 

To summarize, we have shown that the experiments 
involving self-induced transparency are a good alterna- 
tive to the EIT experiments in order to observe optical 
precursors well ahead of the main field, both having in- 
tensities comparable to that of the step-modulated in- 
cident wave. By using a plane-wave model and the rate 
equations approximation, we have established an analyti- 
cal expression for the time-delay of the main-field arrival, 
which generalizes that previously obtained in the infinite 
saturation limit, and we have shown that this expression 
provides a good estimate of the real time-delay as long 
as precursors and main field are well separated and of 
significant amplitude. 
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